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Based on hemagglutinin (HA) and neuraminidase (NA), influenza A virus is divided into 18 dif- 
ferent HA(H1 to H18) and 11 NA types (N1 to N11), opening the possibility for reassortment be- 
tween the HA and NA genes to generate new HxNy subtypes (where x could be any HA and y 
is any NA, possibly). In recent four years, since 2010, highly pathogenic avian influenza (HPAI) 
viruses of H5N1 subtype (HPAI A/H5N1) have become highly enzootic and dynamically evolved 
to form multiple H5 HA clades, particularly in China, Vietnam, Indonesia, Egypt, Cambodia, 
and Bangladesh. So far, after more than 10 years emerged in Vietnam (since late 2003), HPAI 
A/H5N1 is still posing a potential risk of causing outbreaks in poultry, with high frequency of 
annual endemics. Intragenic variation (referred to as antigenic driffi in HA (e.g., H5) has given 
rise to form numerous clades, typically marking the major timelines of the evolutionary status 
and vaccine application in each period. The dominance of genetically and antigenically diver- 
sified clade 2.3.2.1 (of subgroups a, b, c), clade 1 .1 (1 .1 .1/1 .1 .2) and re-emergence of clade 7.1/7.2 
at present, has urged Vietnam to the need for dynamically applied antigenicity-matching vac- 
cines, i.e., the plan of importing Re-6 vaccine for use in 2014, in parallel use of Re-1/Re-5 since 
2006. In this review, we summarize evolutionary features of HPAI A/H5N1 viruses and clade 
formation during recent 10 years (2004-2014). Dynamic of vaccine implementation in Vienam is 
also remarked. 

Keywords: Orthomyxoviridae, Subtypes, HPAI A/H5N1, Genotypes, Clades, Reassortment, 
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Avian influenza is an acute infectious disease that occurs and spreads very fast in poul- 
try, waterfowls, wild birds, animals, and able to zoonotically be transmitted to humans. 
Avian influenza is listed as one of the most dangerous diseases by World Organisation 
for Animal Health (OIE; http://www.oie.int/eng/normes/mmanual/A_summry.htm) 
that a special attention is needed for control. Influenza A virus belongs to the family of 
Orthomyxoviridae, of which the envelope is embedded by two surface glycoproteins, 
e.g., hemagglutinin (HA) and neuraminidase (NA). The virus has a spherical appear- 
ance with a diameter of 80-120 nm, or polymorphological shape, sometimes is elon- 
gated, like fibers and arranged in a sequence. Virus particles (virions) have a simple 
structure, which includes envelope, capsid and a core containing genetic material (vi- 
ral genome) with molecular weight of about 250 million daltons [1] . Influenza A virus 



http://www.ecevr.org/ 



117 



CLINICAL and EXPERIMENTAL VACCINE RESEARCH 

Thanh Hoa Le et al • Dynamic of H5N1 clades and vaccines in Vietnam 



genomes contain a single-stranded negative RNA ([-]ssRNA) 
of 8 segments (polymerase basic 2 [PB2], polymerase basic 1 
[PB1], polymerase acid [PA], HA, nucleoprotein [NP], NA, 
matrix [M], non-structural polypeptides [NS]) [1,2], encoding 
for 15 proteins in overall, including PB2, PB1, PA, HA, NP, NA, 
Ml, M2, NS1, NS2, and supplementary proteins recently iden- 
tified with little or unknown functions [3], i.e., PB1-F2 [4], PB1- 
N40 [5], PA-X [6], PA-N155 and PA-N182 [7] . Up-to-date, there 
are 18 HA and 1 1 NA types for influenza A viruses to be known, 
based on the composition and activity of the HA (from HI to 
H18) and the NA (from Nl to Nil) genes, especially, by re- 
cent addition of two new HA (ie., H17 and H18) and two new 
NA types (N10 and Nil) detected in bats [8,9]. The first 15 HA 
and 9 NA types are originally found in aquatic fowles which 
served as a source for subsequent reassortment [3]. The nu- 
merous types of HA and NA enable the components going 
through genetic reassortments between HA and NA for gen- 
eration of new subtypes [10-12]. 

Major polypeptides encoded by 8 segments of the genome 
of the influenza A virus are as follows: 

1) The polymerase complex (PB2, PB1, and PA): The poly- 
merase enzymes encoded by segments 1, 2, and 3, respec- 
tively, are responsible for transcription and replication of vi- 
ral RNA [13]. Additionally, the recently identified polypep- 
tides of PB1-F2, PB1-N40, PA-X, PA-N155, and PA-N182 en- 
coded within segments 2 and 3, respectively, have confirmed 
the influenza virus as a top of the most complicated viral agents 
by multi-disciplinary utilization of the genome and the func- 
tional polypeptides [3]. 

2) The surface antigenic proteins (HA and NA): HA assists 
virus in attachment to the cell receptors; and by being cleaved 
by protease, facilitating intrusion of viral genetic material into 
the cell. NA functions in promotion of the assembly and the 
release of the virus progenies from the cell receptor [14]. The 
HA and NA antigenic glycoproteins determine the specificity 
of different subtypes of influenza A viruses and NA is a target 
for antiviral treatment [15]. Moreover, HA and NA also play a 
pivotal role in determining the antigenicity of vaccine pro- 
duction [11,14,16]. 

3) The proteins encoded by internal genes: They include 
NP, M (Ml and M2), and NS (comprising NS1 and NS2): Prod- 
ucts of these genes play important roles in packaging of viral 
genomic RNA during replication, in intracellular transporting 
pathways and in shipping mechanism through the channel 
components [17,18]. Ml protein sustains the virion structure, 
playing multiple role in the replication, assembly, and release 
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of viruses [1]. The viral ribonucleoprotein complexes (vRNPs) 
contain viral RNA, viral polymerases and viral NPs which are 
essentially formed during viral infection and replication [19, 
20]. NS2 product also known as a nuclear export protein or 
NEP is essential for transportation of the vRNPs to cytoplasm 
[!]• 

After more than 10 years of highly pathogenic avian influ- 
enza (HPAI) A/H5N1 emergence in Asia, the reassortments 
of HAs and NAs have always occured to generate new sub- 
types [21,22]. Just in a short time of recent 5 years, more le- 
thal, pandemic candidate and/ or pandemic subtypes of in- 
fluenza A viruses have been reassorted and emerged in wide 
geographical localities [3]. These are represented by, namely, 
the H5N1 variants of highly differentiated clades, the pan- 
demic A/H1N1-2009, the new subtypes of H7N9 and H10N8 
in China that posed a high risk across different nations [23]. 
The rapid emergence of the antigenic drifted HPAI H5N1 and 
the newly reassorted H7N9 and H10N8 viruses raise the ur- 
gent needs to definitely update vaccine development from 
antigenicity-matching strains for actual application [15,24- 
26]. Fortunately, a great technology, reverse genetics has con- 
tributed at a great speed, to solve the vaccine development in 
the shortest time, that new vaccine(s) against most recently 
emerging subtypes, with high immunity and warranted pro- 
tection, including the H5N1 and H7N9 vaccines have been 
generated [25-30]. 

In this review, we focus on the recently-introduced HA evo- 
lution, preperably, on H5 HA of HPAI A/H5N1 strains and the 
evolutionary features of clade formation in Vietnam, by 1) 
summarizing the antigenic drift of H5 HA and the emergence 
of its clade/subclades in the past eleven years, since 2003; 2) 
analyzing phylogenetic relationships of different HA clades 
between Vietnam and neighbouring countries; 3) looking at 
spatial-temporal transmission of clades and dynamic vaccine 
implementation in Vietnam. Literature has been taken for use 
mainly from the international sources but in some cases, from 
some of local publications with actual citations. 

Evolutionary Features for Emergence of 
Genotypes and Clades of H5N1 in Vietnam 

It is noted that the precursor highly pathogenic A/H5N1 avi- 
an virus was first time formed in 1996, as a result of the genetic 
reassortment of the H5 HA and the Nl NA types in wild geese 
in Guangdong, China (A/Goose/Guangdong/l/96[H5Nl]); 
and the reassortants subsequently derived from this ances- 
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tral source, have since then spread over the world causing 
numerous epidemics [3,14,31]. Although the more ancient 
strain of A/H5N1 was, indeed, recorded in 1959, indicated by 
the low pathogenic strain of A/ chicken/Scotland/ 1959(H5N1) 
(GenBank: CY015081, GU052518 for HA genes), nearly 20 
years in the past, particularly in the last 10 years, due to the 
actual emergence of HPAI influenza A viruses, the HPAI A/ 
H5N1 subtypes have caused enormous loss to poultry, wild 
birds and humans [3]. Although in history, the evolutionary 
pressure has posed substantial impact on the dynamics of 
the genetic variations of all the 8 genomic segments in order 
to diversify the pathogenicity/ antigenicity and the host adapt- 
ability of the influenza A viruses, most of which have frequent- 
ly occured in the HA and NA genes, resulting in global impact 
by generation of epizootic outbreaks and pandemics [12,16]. 

Antigenic HA polypeptide (i.e., HA, for H5N1, referred to 
as H5) is encoded by segment 4, functionally characteristic to 
the specific binding process to the surface receptor on the 
cellmembrane [14-16,32]. HA has the ability to create genetic 
mutations in different sites (especially, substituting the ami- 
no acids at the protease cleaving site, at the point between 
HA1 and HA2) or to reassort to generate variants, altering 
surface antigenic properties leading to the correlative chang- 
es in immune response [14,31,33]. NA antigenic polypeptide 
of H5N1 (i.e., Nl) encoded by segment 6, is an enzymatic pro- 
tein, cutting off virus-cell receptor molecules linked to sialic 
acid (N-acetylneuramic acid) to timely release virus proge- 
nies from the infected cells during viral infection [14,34] . 

Since the first time appearance in 1996 in Guangdong (Sou- 
thern China), there have been overall 12 genotypes formed 
from the ancestral A/Gs/Gdl/1996 H5N1 strain, including 
GD, A, B, C, D, E, X(Xo-X3), V, Y, W, Z(Z + ) and G [35,36]. In a 
short time after the first emergence in Guangdong, the origi- 
nal genotypes (e.g., GD, A, C, D, E) of the Guangdong lineage 
have entirely disappeared, giving way to the next genotypes/ 
subgenotypes (V, W, Xi, X2, X3, Y, Z, and Z+) evolutionarily 
emerged [37]. The appearance of the Z genotype with strains 
of higher pathogenicity in Southeast Asian countries, was an 
evidence for the "antigenic drift" in genetic constellation of 
the influenza A/H5N1 genome [35,38]. Genotype Z is charac- 
terized by the deletions in NA and NS segments. However, 
some viruses isolated in Southern China (Guangdong, Guangxi, 
and Hunan provinces) have substantially diversified to form 
further genotypes, i.e., Z, V, W, G; and the genotype G was the 
result of reassortant between W and Z [35,39]. In Vietnam, 
genotypes Z and G are dominant and long last since 2005 [40- 



42]. The H5N1 strains isolated in Vietnam during 2004-2006, 
were of the Guangdong sublineage, belonging to the geno- 
type Z. Recendy, the additional genotype G has emerged dur- 
ing 2007-2008; and after 2010, both of Z and G genotypes have 
coexisted up to date [40,43]. Two major groups of subgeno- 
types have been found to exist in Vietnam, subgenotype N 
(North) predominant in the North and S (South) in the South 
[36,41]. Between 2007 and 2010, it was noted that the circu- 
lating strains were of mixed sublineages, i.e., Guangdong-like, 
Fujian-like and Qinghai-like co-existed [41-43]. 

In the recent 10 years, the H5 antigenic gene has high ten- 
dency to drift, generating complicated intragenic changes. As 
results, the newly emerged variants have enhanced their adap- 
tability to multi-host transmissions, helping the HPAI A/H5N1 
viruses to spread throughout the world; hence the level of 
pathogenicity has been exacerbated [32,34]. Since 1996, dur- 
ing 8 years by 2004, H5N1 has diversified to generate 10 dis- 
tinct clades (0-9) due to the antigenic mutational drift in the 
H5 gene, all of which have always originally first time formed 
in China [3]. The emergence of the HPAI A/H5N1 viruses in 
early 2000s has accelerated the H5 gene undergoing onto dy- 
namic evolution in a high speed, generating numerous new 
clades/subclades of the second-, third-, and forth-order sub- 
clades, mosdy within clades 1 and 2 [3,44]. Not all the above 
mentioned clades have been introduced into Vietnam (Fig. 1) 
[41-46], but largely, the rooting clades have further diversified 
to form new serial-order subclades, heavily affecting the avi- 
an influenza occurences in poultry and wild birds, and com- 
plicating the epidemiologic situation in Vietnam [40,42,44, 
47-49]. 

By 2014, for example, clade 2 has diversified into different 
subclades and expanded over the world. The 2.3.4 clade alone 
evolved to the forth order subclades of 2.3.4.1, 2.3.4.2, 2.3.4.3, 
2.3.4.4, 2.3.4.5, 2.3.4.6 over China [50] of which the subclades 
2.3.4.1, 2.3.4.2, 2.3.4.3 have emerged and existed in Vietnam 
for a period of 4 years, during 2007-2010 (Fig. 1). The newly 
emerged clade 2.3.2 appeared more expandable in Vietnam, 
Southeast Asia, North Asia and many countries in Eastern 
Europe from 2009 to 2014, widely detected in wild birds, then 
poultry in China (Qinghai lake region), Mongolia, Russia, Ja- 
pan, Korea, Laos, Hong Kong, Vietnam, Japan, Bangladesh, 
India, Bulgaria and recendy, in Indonesia [42,51-58]. These 
strains are now classified into forth and fifth-order clades with 
an additional letter of the right-most digit implemented [44] . 

Major clades, such as 2.1.3, 2.2.1, 2.3.2, have also undergone 
substantial evolution in HA gene, resulting in a deep diversi- 
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Fig. 1. Emergence timelines of hemagglutinin-based clades of A/H5N1 in Vietnam (only those clades that emerged and identified in Vietnam 
are presented). Data used for making timelines of H5N1 occurence were collected from the following published sources: Nguyen et al. [41,42], 
Creanga et al. [43], World Health Organization/World Organisation for Animal Health/Food and Agriculture Organization (WH0/0IE/FA0) H5N1 
Evolution Working Group [44], Inui [45]; A/B/C marks for clades were taken from the above listed publications; question mark (?) indicates up- 
coming confirmation of clades in 2014. 



fixation, even within clade 2.3.2.1, there have been novel en- 
demic, dynamically differentiated variants of a, b, c subgroups, 
during 2010-2014 in Vietnam. These 2.3.2.1a/b/c viruses, are 
preferably circulating in waterfowles, sometimes with no symp- 
toms, thus are confusingly affecting the vaccination programs 
as well as complicating the epidemiological situation in a num- 
ber of countries in the region [42,43,58,59]. Although six clades 
have been found among the A/H5N1 strains isolated since 
2003 in Vietnam, including 1, 2, 3, 5, 7, 8, only three (clades 1, 
2 and 7) have rooted to spatially/ temporally last and to mono- 
phyletically diversify into discrete subclades of one-, two-, 
three- or forth-order designations (Fig. 1). These clades of 1, 
2, 7 and their digit number-order subclades seemed to un- 
dergo on in-depth H5 based-drifting for formation of the fur- 
ther fifth-order groups, e.g., in case of clade 2.3.2. la/b/c, offi- 
cially nomenclatured by adding alphabetical implementation 
[44]. The derivative strains of these clades are firmly now cir- 
culating in multiple species of gallinaceous poultry through- 
out Vietnam [42,43,48,60]. 

Diversification of the Predominant Clades (1, 
1.1,2.3.4,2.3.2.1 and 7,7.1,7.2) in Vietnam 

Approaching 18 years from the time of first generation, the 
Gs/Guangdong lineage-H5Nl avian influenza virus reassor- 



tants have proved to be the most persistent subtypes express- 
ing the highest pathogenicity in wild birds and domestic poul- 
try [3]. The further descent lineage-viruses spread as largely 
as reaching overall 63 countries with deaths of 250 million of 
gallinaceous poultry and wild birds; and hundreds of human 
victims [25] . Vietnam was among one of the first countries out- 
side of China affected by this HPAI H5N1 subtype and become 
enzootic with a number of persistent outbreaks annually. 

Numerous clades have been generated during last ten years 
indicated by major timelines as below: 

1) During 2003-2005: It was the first time in late 2003, that 
the HPAI A/H5N1 viruses have been introduced into Vietnam; 
and during this period, viruses of clade 1 were predominantly 
circulating in all provinces of the country [40,41] . 

2) In 2006: Vietnam has applied nationwide vaccination 
program (Re-1, clade 0 vaccine), and there was no report of 
any outbreak this year. 

3) Between 2007 and 2009: This period was geographically 
marked by the major presence of clade 2.3.4 of the Fujian sub- 
lineage (and its further subclades 2.3.4.1, 2.3.4.2, and 2.3.4.3) 
in the North; then these (sub)/clades have completely disap- 
peared from the country; no report of clade 2.3.4 was made 
after 2009. Viruses of clade 1 (and its clade 1.1) were domi- 
nating Mekong Delta provinces in the South, and starting 
evolving into clades 1.1.1 and 1.1.2 [41,42]. Re-1 and Re-5 
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vaccines (clade 2.3.4 origin) imported from China were in use 
and a national project of testing NIBRG-14 masterseed (rgA/ 
Vietnam/ 11 94/2004, clade 1, obtained from National Insti- 
tute for Biological Standards and Control [NIBSC, UK]), was 
launched for vaccine production [61,62]. In this period, the 
first detection of a newly designated clade 7 at ports of entry 
in Lang Son province was published [60,63] . 

4) In 2010: This year was the time for the new clades emerged, 
and alternatively, for the appearance of periodically existed 
viruses [42]. The emergence of newly introduced clade 2.3.2 
and further evolution of clade 2.3.2.1 into 2.3.2. la/b was not- 
ed in the North; and the continual evolution of clade 1.1 of 
the Cambodian lineage into 1.1.1 and 1.1.2 occured in the 
South [42,43,47,64]. Re-1 and Re-5 vaccines were in use. 

5) During 2011-2012: The dynamic formation of a, b, and 
particularly, c subgroups of clade 2.3.2.1 was in progress in 
whole country [43,65]; and because of antigenicity-matching 
failure for clade 2.3.2.1b, the Re-l/Re-5 vaccination was ceased 
in some provinces in the North [66] . NIBRG-14 based vaccine 
(namely, NAVET-Vifluvac) was successfully tested, indicating 
reliable efficacy against viruses of the current clades, particu- 
larly when used with immunomodulating adjuvants like P-glu- 
can [61,62,67]. Clades 1.1.1/1.1.2 were in complete dominance 
in the South and Re-1 vaccine was still in use [47,49] . 

6) During 2013-2014: There was a wide circulation of clade 
2.3.2. la,c in whole country and no report of clade 2.3.2.1b 
was made anymore. Concurrently, clades 1.1.1/1.1.2 were 
heading from the Mekong Delta region to Central Vietnam. 
This time, the re-recognition of clade 7.1/7.2 evolving from 
clade 7 was reported. Due to decline in protection and wide 
spread of clade 2.3.2.1c, Re-6 vaccine is planned to import for 
use in 2014 [49]; according to the Vietnam authority, Depart- 
ment of Animal Health, 2014 (http://www.cucthuy.gov.vn/). 

Tracing back the original source of H5N1, it is indicated 
that all viruses of the pre-existing clades 2.3.4, 2.3.2, 1.1, 7 be- 
fore 2010; and the current circulating 2.3.2. la/ (b)/c, clades 
1.1.1/1.1.2 and clades 7.1/7.2 in Vietnam have been introduc- 
ed from Southern provinces of China during 2003-2014 [43]. 
It should also be noted that the very virulent clade 1.1.1/1.1.2 
viruses currently dominating the South of Vietnam have evolv- 
ed from a Cambodian sublineage of clade 1.1 [3]. To examine 
the clustering of viruses in different clades using World Health 
Organization (WHO)/OIE/Food and Agriculture Organiza- 
tion of the United Nations (FAO)-2014-nomenclaturing sys- 
tem, we include those previously listed as references for clad- 
ing source and the current sequences obtained by our study 



from samples collected during 2008-2014. The phylogenetic 
tree resulting from H5 nucleotide data (using 77 sequences 
comprising of representative clades) is shown in Fig. 2. It 
clearly shows that most recent 2009-after collected isolates of 
Vietnam are completely clustered with the H5N1 HPAI strains 
of newly described clade 2.3.2.1 (and of its clade 2.3.2.1a/b/c) 
emerged after 2009, closely linked to A/Hubei/1/2010 (clade 
2.3.2.1a), to A/barn-swallow/Hong Kong/1161/2010 (clade 
2.3.2.1b) and A/Hong Kong/6841/2010 (clade 2.3.2.1c), re- 
spectively. 

Strains forming clade 2.3.2.1a and clade 2.3.2.1b, compris- 
ed HA sequences (majority are of our study, be included in 
analysis) which were closely linked to A/Hubei/1/2010 and 
A/barn-swallow/Hong Kong/1161/2010, respectively. Clade 
2.3.2.1c viruses, all clustered together, around the reference 
A/Hong Kong/6841/2010 strain, are of some typical charac- 
teristics of taxonomic grouping. It is indicated that the old 
strains, originally isolated in 2009 at Qinghai lake (China), or 
those before 2009 including the ancestral A/BHG/Qinghai/ 
1/2009 (accession No. HQ020367) [68] formed a basal branch, 
formerly named as clade 2.3.2, now classified as clade 2.3.2.1c 
(Fig. 2). The tree shown in Fig. 2 supported a precise topology 
of clade 2.3.2.1c strains, with the Vietnamese (of those isolat- 
ed in 2014 of our study as well) and Indonesian viruses [58] in 
two separate groups. Distinct from clade 1 they have evolved 
from, all HA sequences of clades 1 . 1 . 1 / 1 . 1 .2 performed a mono- 
phyletically close group, regardless of where they were identi- 
fied from isolates, in Vietnam or Cambodia (Fig. 2). 

Antigenic Characteristics of A/H5N1 Viruses 
in Vietnam 

Since their re-emergence in 2003, HPAI A/H5N1 viruses have 
become enzootic in some countries and continue causing 
outbreaks in poultry as well as sporadic human infections. 
The A/H5N1 viruses have diversified both genetically and 
antigenically leading to the need for multiple vaccine devel- 
opment. Evolution of the HA gene of influenza A/H5N1 vi- 
ruses emerged in Vietnam has marked major antigenic chan- 
ges, affecting immune response in the vaccinated poultry 
[42,47]. 

Amino acid analysis of H5 HA of the emerging strains iso- 
lated during 2010-2013 in Vietnam has revealed that the poly- 
basic cleavage site (amino acid 338-345, H5 numbering) of 
these viruses, was similar to the H5 HAs belonging to viruses 
of clade 2.3.2.1 and its subgroups (2.3.2. la,b,c) recently re- 
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HA 
(NJ) 




2.3.2.1a 



A-Dk-VN-OIE-2533-2011-AB7O0635 
A-Dk-VN-QT11c4-2012 a > 
A-Dk-VN-QT11c3-2012 a > 

I A-Md-VN-LBM260-2012-AB787209 

A-Dk-VN-OIE-1287-2012-AB735615 
A-Dk-VN-LBM132-2012-AB742263 
A-Dk-VN-NA5c1-2012 a > 
- A-Md-VN-LBM113-2012-AB742255 
A-Ck-VN-KH23-201 3 a| 
A-Ck-VN-TH2c1-2012 a > 
Dk-VN-QT 10-201 2 a} 
A-Ck-VN-TH2c2-2012 B > 
- A-Md-VN-LBM1 4-201 1 -AB741 550 
A-Dk-VN-QT1 604-20 11 a ' 
A-Dk-VN-QT1603C6-2011 a > 
Dk-VN-QT1603C9-2011 a > 
A-Dk-VN-QT801-2011-JN986881 
A-Dk-VN-QT802-2011-JN986882 
A-H ubei-1 -201 0-C Y098758{2 .3.2.1a) 
A-Dk-VN-HT7-2012 a > 
A-Md-VN-LBM57-2011-AB741558 
2 j A-BHG-Mongolia-X53-2009-HM006730 
P- A-BHG-Qinghai-1-2009-HQ020367 
-A-Dk-Laos-469-2010-CY098344 

A-HK-6841-2010-HQ636461(2.3.2.1c) 

50| A-Dk-Sukoharjo-BBVW-1428-9-2012-KC41727 
A-Md-Tegal-B B VW- 1 732-1 1 -201 2-KC41 7275 
A-Dk-Sleman-BBVW-1463-10-2012-KC417273 
A-Dk-Bantul-BBVW-1443-9-2012-KC417272 
A-Dk-VN-HT5-2014 a ' 
^■A-Dk-VN-HT2-2014 a| 
A-Dk-VN-HT4-2014 a > 
A-Md-QuangNinh-6-2013-AB824263 
A-Ck-VN-OIE-2215-2012-AB766239 
A-Dk-VN-OIE-2202-2012-AB769252 
A-Md-VN-LBM227-2012-AB786676 
-A-Dk-VN-KH24-2013 a > 
A-Mdk-VN-LBM436-2013-AB849460 
OOj- A-Md-VN.18153.2009.JN055369(2.3.2.1) 
A-Ck-VN-20346-2009-JN055371 

A-barn-swallow-HK-1161-2010-KC357320(2.3.2.1b) 
A-Md-VN-LBM66-2011-AB741566 
A-Dk-Hunan-S4234-2011-CY146708 
A-Dk-Zhejiartg-213-2011-JN646713 




2.3.2.1c 



2.3.2.1 



2.3.2.1b 




A-Ck-Bangladesh-11rs1984-30-2011-JN795924| 2 3 4 2 
A-VN-HN31432M-2008-HM114617{2.3.4.2) I 
A-Ck-VN-HaTay-44-2007-JX420174 |2 3 4 1 

A-Dk-VN-NCVD-46-2007-CY030548{2.3.4.1)l 
A-Ck-VN-NCVD-40-2007-CY030505 I 
A-Anhui-1-2005-HM172104 2.3.4 
A-Guangxi-1-2005-DQ371930(2.3.4) ' 

A-VN-UT313941-2008-HM1 14593 
A-VN-HN31242-2007-EU294370 

A-Dk-VN -N A72-2007-JX02 1 305(2 .3 
A-Dk-VN-NA114-2007-JX021306 

A-Ck-VN-NCVD-016-2008-FJ842476(7 .1 } I 7 . 1 

A-Ck-VN-NCDV03-2008-FJ842481 



2.3.4.3 




A-Ck-VN-NCVD-093-2008-FJ842480(7.2) 

A-HK 156-1 99 7-AF046088 | Q 



-Gs-GD1-1996-AF148678(0) 

- A-Ck-VN-HD1-2004-EF057807 

- A-VN-1194-2004-EF541402(1) 

-A-Md-VN-GL-2004-EF051515 
-A-Dk-VN-AG-2005-EF051514 

- A-Ck-VN-DT382-2008 « 
- A-Ck-VN-TV98-2008 a > 
- A-Ck-VN-KG88-2008 a > 

- A-Dk-VN-OIE-0062-2012-AB727935(1 .1 .1 ) 
A-Cambodia-S1 21 1 394-2008-HQ200596 
A-Dk-VN-ST0970-2009-JX021303 

A-Dk-VN-OIE-3002-201 1-AB715369 
A-Cambodia-V0203306-2011-JN588805 
A-Dk-VN-OIE-3509-201 1 -AB71 6329 
A-Cambodia-V0219301-2011-JN588806 

A-Md-VN-OIE-0043-2012-AB728619(1.1.2) 
A-Ck-VN-ST-1 -201 2-AB823747 
99 l A-Ck-VN-ST-5-2012-AB818505 



7.2 



1.1.1 



1.1.2 



Fig. 2. Phylogenetic tree based on the hemagglutinin H5 sequences 
(1,588 nucleotides) of representative A/H5N1 viruses including those 
of predominant clades in Vietnam since 2004. Topology was construct- 
ed by MEGA6.06 using the neighbor-joining analysis of the boostrap 
method with 1,000 replications (values <40 were not shown) [72]. 
Clades were nomenclatured based on the criteria and the reference 
strains determined in World Health Organization (WH0)/World Organ- 
isation for Animal Health (0IE)/Food and Agriculture Organization of 
the United Nations (FAO)-2014 system [44]. Major clades of H5N1 are 
shown as groups of monophyletic branches. Regarding to clade 2.3.2.1, 
viruses are divided into clade 2.3.2.1a (A/Hubei/1 /201 0-like); 2.3.2.1b 
(A/barn-swallow/Hong Kong/1 161/201 0-like); and 2.3.2.1c (A/Hong 
Kong/6841/201 0-like). Note: Clades/subclades for prototype reference 
strains are bolded and parenthesized for indication. VN, Vietnam; HK, 
Hong Kong; KR, South Korea; KH, Cambodia. 'Viruses from our own 
study, marked at end of each strain. Scale bar at bottom indicates the 
number of nucleotide substitutions per site; the corresponding Gen- 
Bank accession number are denoted at end of each sequence where 
applicable. 



ported in China, South Korea, Japan, Mongolia, Nepal, Laos, 
Russia, Bangladesh and Indonesia [3,58,69,70]; and similar to 
those of clade 2.3.2 origin, but different from clade 2.3.4 and 
clade 1 (and 1.1) (Table 1) [3,42,43]. The strains found in South 
of Vietnam have evolved from clade 1, resulting in subclades 
1.1.1 and 1.1.2, similar to that of the Cambodian clade 1.1 [3, 
44,70] . However, amino acids of receptor binding site of strains 
(QSG, position 238-240) in all clades in Vietnam and else- 
where were conserved (Table 1), keeping the ability of bind- 
ing specificity stable throughout the evolution of HA [14,59]. 

Analysis of amino acid identity at the antigenic sites (Table 
1) revealed that all strains of clade 2.3.2. la/b/c have amino 
acid K at antigenic position 69; ANPA at 99-102 (for some 
strains, N was replaced by K or S); and a highly conserved se- 
quence, DHEASL, at 140-145 (antigenic site 1) except for the 
Vietnamese clade 2.3.2.1c strains of 2014 who possess NHEVSL. 
A/Hubei/1/2010 and A/Barn-sw/HK/1161(2010), the proto- 
type vaccine strains [30,44], share close identity to the viruses 
within clade 2.3.2.1 at most of major sites, except at antigenic 
site 2, where strains of clades 2.3.2. lb,c were noticed making 
the sequence PYQGNS to become SYQGNS (Table 1). Wheth- 
er these amino acid substitutions, alternatively, at the anti- 
genic site 2, i.e., S152P and K156S/N, between strains of clades 
2.3.2. lb,c and 2.3.2.1a, have caused antigenicity-mismatch- 
ing immune response or not, still a question to be clarified. 

Glycosylation at amino acid 170 in HAi (position 170-172) 
is completely lost by changing from NST (as seen in clade 1 
and 1.1 strains) and NNT (clade 2.3.4) to amino acids DNA in 
all strains of clades 2.3.2 origin and current 2.3.2.1(a,b,c) of 
Vietnam, China, Bangladesh, Indonesia and other countries 
[43,59,71], including those strains collected in 2014 of our 
current study (Table 1). 

Spatial-Temporal Transmission of H5N1 
Clades and Dynamic of Vaccine Application 
in Vietnam 

Clade 2.3.2.1 and clade 1.1 (now 1.1.1/1.1.2) have revealed a 
great diversification due to their complex evolution of HA 
and endemically spatial transmission after 2010 in Vietnam. 
Clades 7 and 7.1/7.2 were temporarily detected before and 
now re-emerged in some places. 

These evolutions and transmissions are characterized by 
the next descriptions, as follows: 

1) Viruses of clade 2.3.2.1a (phylogenetically designated as 
A/Hubei/l/2010[H5Nl]-like), emerged first time in mid-2010 
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Table 1. Amino acid residues at sites of protease cleavage, receptor binding, epitopes of antigenic sites and a site of variable glycosylation of 
the H5 hemagglutinin polypeptide 



Virus strains 


Country/ 


Clade 


Cleavage site 


Receptor 
binding site 






Antigenic sites 




Glycosylation 
change 




Territory 




338-345/346 


238-240 


69 


99-102 


140-145 Site 1 


152-157 Site 2 


170-172 


A /ni/A/M /UTR /9m yl\a) 
A/UK/VINI/H 1 0 (ZU 1 4) 1 


Vietnam 


L.o.l. IC 


nprpppi/ pp 
Unbnnnlvnb 


ncp 
Uob 


y 
h 


A MDA 

AINIrA 


Mur\/ci 
INIblbVoL 


cvnPMC 
oYUblNIo 


nMA 
UINIA 


A /PI/ A/M /^U99 /9H1 Q\a) 

a/lk/vim/khz J (zu i j) ' 


Vietnam 


L.O.I. IC 


nprpppi/ pp 
Unbnnnlvnb 


ncp 
Uob 


Y 
l\ 


A MDA 

AINIrA 


nur aci 
UblbAoL 


cvnPMC 
oYUblNIo 


nMA 
UINIA 


a /ni/ finnii /1 "707 11 nnni 

a/uk/ legal/ 1 /z/- 1 1 (zuiz) 


Indonesia 


Z.J.Z. IC 


nprpppi/ pp 
Unbnnnlvnb 


ncp 
Uob 




A MDA 

AINIrA 


nur aci 
UblbAoL 


cvnPMC 
oYUblNIo 


nMA 
UINIA 


a /pur /ninn!-ni /1 /onno\ 
A/bnb/uingnai/ 1 (zuua) 


China 


9 9 9 1 o 

Z.J.Z. IC 


nprpppi/ pp 
Unbnnnlvnb 


ncp 
Uob 


Y 
l\ 


A MDA 

AINIrA 


nur aci 
UblbAoL 


DVPPMC 

rYUblNIo 


nMA 
UINIA 


A/UI//CQ/11 /omn\b) 
A/nl\/oo4 1 (ZU IU) ' 


Hong Kong 


9 9 9 1 o 

Z.J.Z. IC 


nprpppi/ pp 
Unbnnnlvnb 


ncp 
Uob 


Y 


A MDA 

AINIrA 


nur aci 
UblbAoL 


cvnPMC 
o i UbINIo 


nMA 
UINIA 


A /l\/lrll/A/M /I pmgg /9m 1 \ 
A/MuK/ VIM/LDlVIbb (ZU I I ) 


Vietnam 


9 Q 9 1 h 

Z.J.Z. I D 


nitppppk'pp 
UlbnnnnKnb 


ncp 
Uob 


Y 

Is 


A MDA 

AINIrA 


nur aci 
UblbAoL 


cvppmc 
oYUblNIo 


nMA 
UINIA 


a /Dim niAf/u^/1 1 ci /onm\b) 
A/barn-sw/HK/ 1 loi (ZUlU)' 


Hong Kong 


9 9 9 1 h 

L.O.L. I D 


nmpppp^pp 
UlbrinnnKhb 


ncp 
Uob 


Y 


A MDA 

AINIrA 


nur aci 
UblbAoL 


cvnPMC 
oYUblNIo 


nMA 
UINIA 


a /ni/A/M /i pmi 9G /9m 9\ 
A/UK/VINI/LdIVI I db (ZU IZ) 


Vietnam 


9 9 9 1 t 

L.O.L. la 


nprpppi/ pp 
Unbnnnlvnb 


ncp 
Uob 


Y 
l\ 


A MDA 

AINIrA 


nur aci 
UblbAoL 


DVPPMC 

rYUblNIo 


nMA 
UINIA 


A /PUA/M /TUOoi /onnte) 
A/LK/VINI/ 1 HZC I (ZU I Z) ' 


Vietnam 


9 9 9 1 t 

z.j.z. la 


nprpppi/ pp 
Unbnnnlvnb 


ncp 
Uob 


Y 


A MDA 

AINIrA 


nur aci 
UblbAoL 


DVPPMC 

rYUblNIo 


nMA 
UINIA 


a /ni/ A /m /PTom /omi\al 
A/UK/VIM/UIdU I (ZUI I) 


Vietnam 


9 9 9 1^ 

z.j.z. la 


PDrDPDI/ DP 

Unbnnnlvnb 


ncp 
Uob 


Y 


A MDA 

AINIrA 


nur a ci 
UnbAoL 


Dvnpoc 
rYUbbo 


nM a 
UINIA 


A/HUDei/l (ZUlU) 


China 


9 9 9 1/1 

L.o.L. la 


PDCPPDI/ DP 

Unbnnnlvnb 


ncp 
Uob 


Y 

K 


A MDA 

AINIrA 


nur a ci 
UnbAoL 


Dvnpi/c 
rYUbKo 


nM a 
UINIA 


A/ft/liinL/L^D/l/in/cic/onii\ 
A/MalUK/Rn/K I U/b I b (ZU I I ) 


South Korea 


9 9 9 1 

L.o.L. I 


PDCPPDI/ DP 

Unbnnnlvnb 


ncp 
Uob 


Y 


A l/DA 

AI\rA 


rue A CI 

bnbAoL 


DVPPMC 

rYUblNIo 


nM a 
UINIA 


A /Df/ilo/^r, /Tnoliir^i /1 C /9H1 1 \ 

A/naicon/ locnigi/ lb (ZU I 1) 


Japan 


9 9 9 1 

Z.J.Z. I 


nprpppi/ pp 
Unbnnnlvnb 


ncp 
Uob 


Y 
K 


AI/DA 

AI\rA 


nur aci 
UblbAoL 


cvppmc 
oYUblNIo 


nMA 
UINIA 


A /P/iPr-/^k/i /nin^knl /1 /OHnflV 

A/bcbrene/Uingnai/ 1 (ZUUa) 


China 


9 9 9 1 

L.o.L. I 


nDrppDi/ dp 
Unbnnnlvnb 


ncp 
Uob 


Y 


A MDA 

AINIrA 


nur a ci 
UnbAoL 


DVPPMC 

rYUblNIo 


nM a 
UINIA 


A /DUP /ft Annnnlin /VC9 /9f"inn\ 

A/b>nb/lVlongolia/Abj (zUUa) 


Mongolia 


9 9 9 1 

L.o.L. I 


nDrppDi/ dp 
Unbnnnlvnb 


ncp 
Uob 


Y 

K 


A MDA 

AINIrA 


nur a ci 
UnbAoL 


CVPPMC 

oYUblNIo 


nM a 
UINIA 


A/PL/MD/01 /omn\ 

A/L.K/IMr/0 1 (ZU I U) 


Nepal 


9 9 9 1 

L.o.L. I 


nDrppDi/ dp 
Unbnnnlvnb 


ncp 
Uob 


Y 


A MDA 

AINIrA 


nur a ci 
UnbAoL 


cvnPMC 
MUblNlo 


nM a 
UINIA 


A/nu/i a //icq /omn\ 
A/UK/LA/4by (ZU IU) 


Laos 


9 9 9 1 

Z.J.Z. I 


nprpppi/ pp 
Unbnnnlvnb 


ncp 
Uob 


Y 
K 


A MDA 

AINIrA 


nur aci 
UblbAoL 


CVPPMC 

oYUblNIo 


nMA 
UINIA 


A ffmkn /Ti n ,/i /9 /OOnflV 

A/brebe/ lyva/j (zuuy) 


Russia 


9 9 9 1 

L.o.L. I 


nDrppDi/ dp 
Unbnnnlvnb 


ncp 
Uob 


Y 


A l/DA 

AI\rA 


nur a ci 
UnbAoL 


cvnPMC 
oYUblNIo 


nM a 
UINIA 


A/Littlebgret/rlK/oobj (zuu/) 


Hong Kong 


9 9 9 

L.o.L 


nDrppDi/ dp 
Unbnnnlvnb 


ncp 
Uob 


Y 
h 


A MDA 

AINIrA 


nur a ci 
UnbAoL 


DVPPMC 

rYUblNIo 


nM a 
UINIA 


A /PnnUnrnn /UY /C^Q /9nOQ\ 

A/bcnbieron/niybba (ZUUo) 


Hong Kong 


9 9 9 

L.o.L 


nprpppi/ pp 
Unbnnnlvnb 


ncp 
Uob 


Y 
l\ 


A MDA 

AINIrA 


nur aci 
UblbAoL 


DVPPMC 

rYUblNIo 


nMA 
UINIA 


A /Pi mnnui /1 /OnflQ\ 

A/buangxi/ 1 (ZUUo) 


China 


9 9/19 

Z.J. 4. J 


i prpppi/ pp 
Lnbnhnlvnb 


ncp 
Uob 


Y 


A MDA 

AINIrA 


nur acc 
UblbAoo 


DvnPCD 

rYUbor 


MMT 

INI INI I 


A A/M /UM91 9/19 /onn7\ 
A/VINI/HIN|jlZ4Z (ZUU/) 


Vietnam 


9 9/19 

Z.J.4.J 


i prpppi/ pp 
Lnbnhnlvnb 


ncp 
Uob 


Y 
l\ 


A MDA 

AINIrA 


nur acc 
UblbAoo 


DVPP\/D 

rYUbVr 


MMT 

INI INI I 


A/Ck/VN/NCVD-093 (2008) 


Vietnam 


7.2 


QIEGRRRKRG 


QSG 


K 


ASPA 


NHETSL 


SYLENP 


NNT 


A/Ck/VN/NCVD-01 6(2008) 


Vietnam 


7.1 


QREGGRRRKG 


QSG 


K 


ASPA 


NHEASS 


SYLGEP 


NNT 


A/MdkA/N/OIE-0043(2012) 


Vietnam 


1.1.2 


QREERRKKRG 


QSG 


R 


ANPV 


SHEASL 


PYQGQS 


NST 


A/Md/VN/OIE-559 (201 1 ) 


Vietnam 


1.1.1 


QREGRRKKRG 


QSG 


R 


ANPV 


SHEASL 


PYQGNS 


NST 


A/KH/R0405050 (2007) 


Cambodia 


1.1 


QREGRRKKRG 


QSG 


R 


ANPV 


SHEASL 


PYQGNS 


NST 


A/VN/1 194 (2004) 


Vietnam 


1 


QRERRRKKRG 


QSG 


K 


ANPV 


SHEASL 


PYQGKS 


NST 


A/VN/1 203 (2004) 


Vietnam 


1 


QRERRRKKRG 


QSG 


K 


ANPV 


SHEASL 


PYQGKS 


NST 


A/Gs/Guangdong/1 (1996) 


China 


0 


QRERRRKKRG 


QSG 


R 


ASPA 


NHDASS 


PYHGRS 


NSA 



Bold letter of amino acid indicates a change in comparison to other strains of same clade. 

Dk, duck; Ck, chicken; BHG, black headed gull; HK, Hong Kong; Mdk, muscovy duck; VN, Vietnam; sw, swallow; MalDk, mallard duck; KR, South Korea; KH, Cambodia; 
Gs, goose. 

a, Viruses from our own study. 

b| Prototype reference strains for clade 2.3.2.1 a/b/c. 



in some provinces of the North of Vietnam. During the past 
four years since then, HA 2.3.2.1a viruses spatially spread over 
the coastal provinces along Central Vietnam, reaching Me- 
kong Delta Basin, and currently (early 2014) predominantly 
are present in whole country [42,44] . For vaccination of poul- 
try, Re-5, an imported vaccine from China for use in Vietnam, 
also known as Harbin Re-5 (reverse genetics-based A/ duck/ 
Anhui/1/2006 of HA clade 2.3.4 antigens), produced by the 
Harbin Veterinary Research Institute, Heilongjiang province, 



China, is still in nationwide use, showing reasonable efficacy 
for protection against clade 2.3.2.1a viruses [49]. The fact is 
that, circulating viruses of this 2.3.2.1a clade have not caused 
great loss to poultry and humans in Vietnam. Both Re-5 and 
Re-1 vaccines (Re-1, implemented since 2006; of rgA/Goose/ 
Guangdong/1996 of HA clade 0, generated by reverse genet- 
ics) conferred a high level of protection (over 90%) in chick- 
ens against clade 2.3.2.1a by challenge test in laboratory; and 
about 70% seropositivity (by HA inhibition test) in immuniz- 
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ed poulty in the field. However, protection level of Re-l/Re-5 
vaccines, by challenging with strains of novel clade 2.3.2.1 (a, 
b, c), was much lower in ducks and muscovy ducks [47,49,65, 
67,69,73]. 

2) Viruses of clade 2.3.2.1b (A/barn-swallow/HK/1161/2010 
[H5Nl]-like), were first time identified in early 2011 in some 
localities in the North Vietnam such as Quang Ninh, Bac Ninh, 
Bac Kan, Thai Nguyen, Phu Tho, Nam Dinh, Thai Binh, Son 
La, Ninh Binh, Nghe An [42,43,64,66]. Up-to-date, viruses of 
this clade 2.3.2.1b circulated restrictly in China, Hong Kong 
SAR and Vietnam [44]. In Vietnam, clade 2.3.2.1b viruses seem- 
ed to limit their spatial and temporal transmission in the North 
of Vietnam. Nghe An province was the last geographic locali- 
ty where the latest clade 2.3.2.1b was detected in 2012. In 2013, 
no further report has been made about the presence of HA 
clade 2.3.2.1b viruses, assuming that they have completely 
disappeared from the circulating A/H5N1 population in the 
country [44]. 

3) Viruses of clade 2.3.2.1c (the reference sttain therein is 
A/Hong Kong/6841/2010[H5N1]), was first time identified in 
2012, from outbreaks of chickens and ducks in several bound- 
ary northern provinces, then spread southbound to Quang 
Ngai, Khanh Hoa provinces of Central Vietnam. These virus- 
es are now shifting towards the provinces in Mekong Delta of 
Southern Vietnam, and predominantly present in the whole 
country [43,66,67,70]. Thus, clade 2.3.2.1c is of a great chal- 
lenge to HPAI A/H5N1 control in Vietnam, since the current 
in-use-vaccines showed rapid decline of protection against 
these viruses; and of course, presenting a high risk to public 
health [49]. Regarding to nationwide immunization of poul- 
try, the current-in use Re-1 and Re-5 vaccines do not have 
sufficient efficacy to prevent the 2.3.2.1c clade viruses [49], 
requesting a new vaccine with higher matching of antigenici- 
ty to be implemented [43]. The Re-6 vaccine (rgA/dk/Guang- 
dong/S1322/2010[H5Nl] of clade 2.3.2), developed by the 
Harbin Veterinary Research Institute (China), has been im- 
ported for use in 2014, targeting the prevailing clade 2.3.2.1 in 
Vietnam. 

4) Currentiy, clade 2.3.2.1 tends to spread in spatial distance 
from North to South of Vietnam, warning the vertical trans- 
mission along the country. Concurrently, strains of more vir- 
ulent clade 1.1, now evolved into 1.1.1 and 1.1.2, are predom- 
inantly present in Mekong Delta provinces, which have been 
phylogenetically revealed to share close evolution of the Cam- 
bodian clade 1 progenitors [42,44,73]. The presence of clade 
1.1 and its subclades 1.1.1/1.1.2 in several provinces of Cen- 
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tral Vietnam in 2013, indicated that now these clades are head- 
ing towards North Vietnam from Mekong Delta Basin. The 
overlapping region for spatial transmission of clades 1.1/1.1.1/ 
1.1.2 and 2.3.2.1a/b/c is lower Central Vietnam provinces; 
and interestingly, contrary orientation of HPAI A/H5N1 in- 
fections is likely of a consequence of free poultry trade due to 
North-South/South-North transportation, between provinces 
within Vietnam and between Vietnam and Cambodia [43,47]. 

5) Clade 7, 7.1/7.2 viruses (A/chicken/Shanxi/10/2006-like), 
representing an interesting evolution of clades in Vietnam, 
were on site detected in chickens seized at ports of entry in 
2008 [60] and subsequently, viruses of these subclades (7.1/7.2) 
were identified in some places of Northern Vietnam [42,43,63]. 
The early introduced A/chicken/Vietnam/NCVD-016/2008 
(No. FJ842476) was considered as an ancestral entry strain for 
diversification of the recently recognized clade 7.1 and 7.2 vi- 
ruses. 

Conclusion 

From 2009 to date, in recent 5 years, a number of new anti- 
genic subgroups of H5N1 have been formed mainly from two 
sources in Vietnam: newly introduced clades (2.3.2.1) from 
Southern region of China, and clade 1.1 localized in the Me- 
kong Delta evolved from Cambodian lineages, heading to 
move in the opposite direction between the North and South 
Vietnam. Meanwhile, the in-depth diversified clade 2.3.2.1 
viruses (with a, b and c monophyletic groups) have compli- 
cated the epidemiologic and vaccination status in the whole 
country. Because of failure in immune response in poultry 
vaccinated with Re-1 vaccine (clade 0) and Re-5 (clade 2.3.4) 
to viruses of the currently circulating clade 2.3.2. lb/c, the 
Vietnamese government has announced a temporary sus- 
pension of the vaccination program with these vaccines in 
some of the northern provinces. Significant amino acid sub- 
stitutions in the strains of clade 2.3.2.1c emerged since 2011, 
were the main cause of the protective decline of the available 
vaccines currently used nationwide. The biased efficacy of 
the Re-1 and Re-5 imported vaccines has forced Vietnam to 
import Re-6 vaccine (clade 2.3.2) currently being used in 2014. 
The transmission tendency of two predominant clades, i.e., 
2.3.2.1 (a,b,c) and 1.1/1.1.1/1.1.2, needs special works for ac- 
tive management of the endemics under control and con- 
taining enforcement of the virus spread. 
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